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Abstract: We derive the post-Newtonian dynamics for a spinning body with Corinaldesi- 
Papapetrou spin supplementary condition in Kerr spacetime. Both the equations of motion 
for the center-of-mass of body and the spin evolution are obtained. For the non-relativistic 
case, our calculations show that the magnitude of spin measured in the rest frame of the 
body’s center-of-mass does not change with time, though the center-of-mass does not move 
along the geodesic. Moreover, we find that the effects of the spin-orbit and spin-spin 
couplings will be suppressed by the Lorentz factor when the body has a relativistic velocity. 
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1 Introduction 

The dynamics of the point particle in Kerr field has been presented in many classical 
textbooks [1-4], More further studies can be found in recent papers [5-7]. The dynamics of 
the extended body is much more complex, because the internal structure of the body usually 
subjects to the gravitational field and makes the body follow a nongeodesic trajectory. 

The first attempt to explore the dynamics of a spinning body in the curved spacetime 
is performed by Mathisson [8], who shows there is an interaction term coupling the body’s 
spin tensor to the Riemann curvature tensor in the equations of motion. Papapetrou [9] 
derives the equations of motion for the extended body with respect to the proper time, and 
demonstrates that the body moves along the geodesic in the single-pole approximation and 
the spin effect appearing in the pole-dipole approximation makes the worldline of the body 
deviate from the geodesic. The extended body under the pole-dipole approximation is also 
called as the spinning particle, and the general equations of motion for the spinning particle 
are commonly referred to as Mathisson-Papapetrou (MP) equations. 

MP equations are not a closed system, and the spin supplementary condition (SSC) 
needs to be imposed to fix the representative worldline of the extended body. The commonly- 
used SSCs are Corinaldesi-Papapetrou (CP) [10], Pirani (P) [11] and Tulczyjew (T) [12] 
conditions. CP condition chooses the center-of-mass of the body in the rest frame of the 
gravitational source as the representative point. P and T conditions choose the center- 
of-mass of the body in its own rest frame and in the zero 3-momentum frame as the 
representative point, respectively. Although different SSCs lead to different solutions of 
MP equations for the same physical motion, the choice of SSCs still remains a not well 
understood problem [13, 14]. 

The motion of the spinning particle in the field of black hole has been studied by 
many researchers [15-33]. Specifically, Schiff derives the equations for the spin evolution 
in Schwarzschild field under CP and P conditions respectively [15]. Barker and O’Connell 
study the nongeodesic motion of the spinning particle in Schwarzschild field under CP and 
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P conditions, and demonstrate that the equations of motion under these two conditions 
are consistent with each other by shifting in the center of mass [16]. Bini, Gemelli, and 
RufHni invesigate the behaviour of charged spinning particle moving along circular orbits 
in the equatorial plane of Reissner-Nordstrom spacetime under CP, P and T conditions 
respectively [18]. Kyrian and Semerak solve MP equations with various SSCs in Kerr field 
numerically, and show that the trajectories of the representative point of the spinning par¬ 
ticle under different SSCs are different [19]. Plyatskon and Fenyk investigate the highly 
relativistic circular orbits of a spinning particle in Schwarzschild field under P and T condi¬ 
tions [22], and the different cases of the spin orientation and the direction of the particle’s 
orbital motion on the highly relativistic circular orbits in the equatorial plane of Kerr source 
under P condition [23]. Recently, Cralla and Wald also derive the gravitational self-force 
of the spinning particle in Schwarzschild-de Sitter spacetime under CP condition [34]. In 
this work we study the dynamics of the spinning particle in Kerr field under CP condition, 
which so far has not been explored. Basing on the harmonic metric for Kerr spacetime [7], 
we derive the post-Newtonian equations of motion for the spinning particle including the 
spin evolution, to the next-to-leading terms of the spin-curvature coupling. 

The content of this paper is arranged as follows. Section 2 briefly introduces MP 
equations, and gives the equation of motion for the spin under the CP condition in the 
post-Newtonian approximation. In Section 3 we derive the equations for the particle’s 
motion and spin evolution in Kerr spacetime, and calculate the spin precession measured 
in the rest frame of the particle for the non-relativistic case. Summary is given in Section 
4. 


2 Mathisson-Papapetrou equations 


A spinning particle deviates from the geodesic motion and moves along an orbit described 
by MP equations of motion [9] 


D / D N 1 


Ds 


+ vf^Up- 


DS^P 

IDs 


DS»P 


— - = 0 , 


Ds 


( 2 . 1 ) 

( 2 . 2 ) 


where = dx^/ds and 5"^ are the particle’s contravariant 4-velocity vector and spin 
tensor, respectively, m denotes the kinematical rest mass of the spinning particle, and 
it reduces to the usual rest mass when the spin vanishes. D/Ds denotes the covariant 
derivative with respect to the particle’s proper time s. Rx„p is the Riemann curvature 
tensor corresponding to the background on which the particle moves. For the spinless 
particle one can easily verify that the above trajectory Eq. (2.1) reduces to the usual 
geodesic equation. Here and in the following, we use the geometrized units {c = G = 1) 

and the metric signature (-h), with Creek indices running from 1-4 and Latin indices 

from 1-3. 


- 2 - 






The kinematical rest mass m of the spinning particle may not be a constant, and its 
evolution is described by [33] 
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where p°‘ is the contravariant 4-momentum vector of the spinning particle 
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Substituting Eq. (2.3) into Eq. (2.1), we have 
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General speaking, the acceleration of the spinning particle is given by 


dv^ 

dt 


1 


,0'i2 


du^ 

ds ds J 


du^\ _ j j 

I ®geo ~k ®spin ) 


where v'^ = dx''/dt = ud ju^ is the velocity of the particle, and 
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in which Ogg^ and a*pj^, respectively, denote the acceleration caused by the curvature of 
spacetime and the coupling of spin-curvature. 

For the evolution of spin, since is an anti-symmetric tensor, it has six independent 
components. However, Eq. (2.2) only has three independent equations, which can be re¬ 
written as [9] 

DS//^ _u^DS^ _'tDDS^ 

Ds Ds u^ Ds ’ 

therefore, some supplementary conditions have to be imposed. In this work we only consider 
the Corinaldesi-Papapetrou condition [10] 


5*^ = 0 , 


(2.9) 


which implies that x°‘ is the center-of-mass of the spinning particle, and thus makes Eq. (2.1) 
describe the trajectory for the particle’s center-of-mass in the curved spacetime [15]. 

In the post-Newtonian approximation, the affine connection can be expanded as 
follows [1] 
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where the symbol denotes the term in of order /f, with v and f representing 
respectively the typical values of velocity and distance in a non-relativistic system. 

Substituting Eqs. (2.9) and (2.10) into Eq. (2.8), keeping the coefficient of on the 
right-hand side of Eq. (2.8) to order v^/f, we have 


Ds 

Ds 


3 2 

Ttj ni+ rU . 

(2.11) 

rt^+rt4)(rfs“-«*s>'=). 

(2.12) 


Eqs. (2.4), (2.11), and (2.12) are the basis for us to obtain the equation of motion for 
the center-of-mass of the spinning particle, as well as the spin evolution. 


3 Equations of motion in Kerr spacetime 


The spacetime for a constantly rotating gravitational source is referred to as Kerr spacetime. 
In the harmonic coordinates, Kerr metric can be written as [7, 35] 
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where M and a denote the mass and the angular momentum per mass for the grav¬ 
itational source, respectively. The angular momentum is assumed to be codirectional 
with the positive direction of z-axis. x = {x, y, z), x ■ dx = xdy + ydy + zdz, and 
{x‘^ + y‘^)/{r‘^ + 0 ^) + z‘^/r"^ = 1. Expanding the coefficients of the metric Eq. (3.1) in powers 
of 1 /r to next-to-next-to leading order, we have 


9i4 = —(1 + </>)Cj ) 

944 = 1 + 2i^ -|- 2</>^ , 

where (j) = —M/r is Newtonian potential, ^ = 2aM{x x e^)/r‘^ denotes the vector potential 
due to the gravitational source’s rotation, with 63 being the unit vector of z-axis. 
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The corresponding affine connection which are needed in the following calculations can 
be written as follows 


^ • Ml. 

2 2 M ^ I 

Ti,=Tt=^x\ r),=-iCij - CjA 
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Following Schiff’s 


M^ / • ^ 2x^x^x^ 
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® • 3M M 

approach [15], we can define the purely spatial vector S as 

S = {SuS2,S3) = iS^^S^\S^^) . 


(3.2) 


(3.3) 


Plugging Eq. (3.2) into Eq. (2.12), and making use of Eq. (3.3), we can obtain 

^ (x ■ S) V — 2 (v ■ S) X + 3(v ■ ®)5l + -S X (V X C) • (3.4) 

at ■'2 

When the gravitational source’s rotation vanishes, Eq. (3.4) reduces to the result of Schiff 
[15]. 

It is worth emphasizing that the spinning particle does not move along the metric 
geodesic, and Eq. (3.4) is different from the spin evolution given by Eq. (9.6.5) in the 
textbook [1], in which the spinning particle is treated as a point particle following the 
geodesic. 

Following the same procedure given by Schiff, we can obtain the spin evolution mea¬ 
sured in the rest frame of the particle as follows (see Appendix): 


with 


dS, 


rest 


dt 


— n X iSrest ) 


^ . 3M, 

n = --v xc-^ivxx) 
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(3.6) 


The first term in the right hand side of Eq. (3.6) is just Lense-Thirring precession, and the 
second one is called as the geodetic precession or Thomas precession caused by gravita¬ 
tion [1]. 

Although the formula for the precession frequency is the same as Eq. (9.6 .12) in the 
textbook [1], both the equations of motion for the spinning particle and the evolution of the 
spin here differ from those in the textbook. These differences are due to the fact that we start 
with the Mathisson-Papapetrou equations, which are basing on the dipole approximation, 
to derive the equations of motion of the spinning particle and the spin evolution, while in 
the textbook the single-pole approximation is used to obtain the corresponding equations. 
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We can see from Eq. (3.5) that S^est precesses at a rate |ri| around the direction of fi. 
This implies that the magnitude of spin measured in the rest frame of the particle does not 
change with time in the post-Newtonian approximation, though the spinning particle does 
not move along the geodesic. 

Now we turn to derive the equation of motion for the center-of-mass of the particle in 
Kerr spacetime. When taking into account of the effects of spin on the particle’s motion, 
we only keep the coefficient of on the right-hand side of Eq. (2.4) to the next-to-leading 
terms. For the non-relativistic case in which the particle is assumed to be bounded by the 
gravitational source, Eq. (2.4) can be approximated as 
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Substituting Eqs. (3.7) and (3.8) into Eq. (2.5) and making the approximation 1/u'^ ~ 
1 — M/r — w^/2, we can obtain the acceleration of the non-relativistic particle as follow 

dv 3M / 2M 1 2\ rr / 2/ r,/ 

— = aiow H- 5(1 - /'ft® ■ ^ X S) + 2{v ■ x){x x S)} 

-T{(» . V)"(S X C) + (« ■ V)(S X V)(» ■ C) + [»■ V(s. (V X C))|k} 

+ ^^\s X C- SxV{x-C)-V[x-{S X C)] + 2lx ■ (S X V)]C 
2mr^ I 

+3(a;-V)(SxC)+3[S-(VxC)]a;}, (3.9) 

where aio^ represents the acceleration of the non-relativistic particle along the geodesic of 
the Kerr spacetime, which can be found in our previous work [7]. We have also verified that 
when the reference worldline of the particle is shifted to x* = x* — S'^^Uj/m, this equation 
confirms the relative one-body equation of motion in a two-body system [36] in the extreme- 
mass-ratio limit (see Eqs. (2.1) and (2.2) with r/ —)• 0 therein), but our result includes 
higher-order effects of the spin. 

For the spinning particle with a relativistic velocity, we have 0{v^) ~ 1. The evolution 
equations of spatial components of can be derived from Eq. (2.12), which reads 

^ ^ \ {x ■ S) V — 2 {v ■ S) X + 3{v ■ x)5l -|- -5 X (V X ^) 

dtr-^ ■'2 

-[u • V(u • C)]S + ^(^; • 5) [(u • V)C + V(u • 0] . 
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The equation of motion of the relativistic particle can be derived from Eq. (2.4), which can 
be approximated as 
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Substituting Eqs. (3.11) and (3.12) into Eq. (2.5), we can obtain the dynamics for the 
relativistic particle with spin as follow 

^ = “high + ^ X -S') + 2(u ■x)(x X S)} 

at 'ymr° 

-^{(^ ■ X 0 + (^ • V)(5 X V)(u • 0 + [^^ • V(S • (V X C))]^ 

-V[5-(VxC)]}, (3.13) 

where Ohigh is the acceleration of the relativistic particle along the geodesic of the Kerr 
spacetime and has been derived in Ref. [7]. 7 = (1 — is Lorentz factor. It is observed 

that the effects of spin-orbit and spin-spin couplings are suppressed by a factor of 7 in the 
relativistic case. 


4 Summary 

The motion of the spinning particle under the CP condition in the Kerr spacetime has been 
investigated. We derive the post-Newtonian equation for the evolution of spin, and the 
acceleration of the spinning particle due to the coupling of spin-curvature to the next-to- 
leading terms. We also calculate the precession of spin measured in the rest frame of the 
particle when the particle is non-relativistic, and hnd that the measured spin’s magnitude 
does not change with time though the spinning particle does not move along the metric 
geodesic. Moreover, it is found that the effects of the spin-orbit and spin-spin couplings in 
the equations of motion are suppressed by the Lorentz factor when the spinning particle 
has a relativistic velocity. 

A The evolution of spin in the rest frame of the particle. 

Let Srest denote the spin with respective to the rest frame of the particle. When the 
particle’s velocity v is non-relativistic, we have [15] 

5rest= (1 + ^) 5 +^ [u25-(u-5)u] , (A.l) 

where v = lul. 
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Differentiating this equation with the time t, we can get 


; dS dv 
dt dt 


dSr-est if dv ^ 

dt 2\_ dt 

2M dr 
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From Eq. (3.4) we have 


Making use of Eq. (A.3) and the following relations 
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dr/dt = {x ■ v)/r, dv/dt ^—Mx/r^, 
dv/dt = {v ■ dv/dt)/v , 

we can obtain the rate of change of Surest to order (r;^/r)S'rest from Eq. (A.2) as follows 

^ ^ ■ 

Plugging Eq. (3.4) into the above equation, we have 

d'S'rest 


dt 


^ \sx{VxC)+‘^^[{x-S)v-{vS)x] 

= Sx xC + ^{vxx) 

, 3M, A ^ 

X C + ^(r? X *) xSrest , 


(A.4) 




where D represents the precession frequency. 
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